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Abstract 

Dirac deformation of Poisson operators of arbitrary rank is considered. 
The question when Dirac reduction does not destroy linear Poisson pen- 
cils is studied. A class of separability preserving Dirac reductions in the 
corresponding quasi-bi-Hamiltonian systems of Benenti type is discussed. 
Two examples of such reductions are given. This paper will appear in J. 
Phys. A: Math. Gen. 
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1 Introduction 

Recently a new (quasi)-bi-Hamiltonian separability theory of Liouville integrable 
finite dimensional systems was constructed PP-QHI based on general properties 
of Poisson pencils on manifolds. A natural further step within this theory is to 
investigate admissible integrable and separable reductions of these integrable / 
separable systems onto appropriate submanifolds. The most natural approach 
seems to be the one based on the Dirac theory of constrained dynamics |16 ) .|17 ) . 
The presented paper contains only some special cases of such reductions, but 
even in these cases the problem is far from being trivial. The difficulties we met 
during our research inclined us to reconsider the Dirac formalism from the point 
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of view of Poisson bivectors rather then from the point of view of constrained 
dynamics. 

The paper is organized as follows. In this introductory part we remind ba- 
sic concepts of Poisson geometry and of (quasi-)bi-Hamiltonian systems. In 
Section 2 we formulate the theory of Dirac reductions of Poisson brackets in 
terms of Poisson bivectors. Our construction is more general than usually met 
in literature since we consider reduction procedure on the whole foliation of 
submanifolds. Here we also explain that our approach indicates that Dirac 
classification of constraints onto those of "first" and of "second class" requires 
further discussion as our further results show. In Section 3 we review the re- 
cent results on separability theory of quasi-bi-Hamiltonian systems of Benenti 
type. In section 4 we perform Dirac reduction of Poisson pencil and corre- 
sponding quasi-bi-Hamiltonian chain of Benenti type onto a particularly chosen 
submanifold. Chosen constraints preserve the Liouville integrability as well as 
the coordinates of separation of considered system. The main obstacle of such 
a choice is that these constraints are nonexpressible in natural (original) coor- 
dinates. Hence, in Section 5 we modify constraints introduced in the previous 
section and obtain an equivalent reduction, that is expressible directly in origi- 
nal coordinates. Finally, in Section 6 we illustrate the results by two nontrivial 
examples of constrained separable dynamics. 

Let us first remind few basic facts from Poisson geometry. Given a manifold 
A4, a Poisson operator tt on is a mapping tt : T*A4 TA4 that is fibre- 
preserving (i.e. tt\t^m ■ T*M T^M for any x £ M) and such that the 
induced bracket on the space C°°{M.) of all smooth real- valued functions on M 



{.,.}, ■.C°^[M)xC^{M)^C°^{M) , {F,GK - (dF,^dG) (1) 



(where (., .) is the dual map between TAi and T*Ai) is skew-symmetric and sat- 
isfies Jacobi identity (the bracket always satisfies the Leibniz rule {F, GH}^ = 
G {F, H}^ + H {F, G}^). Throughout the whole paper the symbol d will denote 
the operator of exterior derivative. The operator tt can always be interpreted 
as a bivector, tt G A-^(A^) and in a given coordinate system (cc^, . . . , x™) on 
we have 



A function G : ^ R is called Casimir function of the Poisson operator tt 
if for arbitrary function F : ^ R we have {F, G}^ = (or, cquivalently, if 
TrdG = 0). A linear combination tt^ — tti — ^ttq (^ G R) of two Poisson operators 
TTo and TTi is called Poisson pencil if the operator tt^ is Poisson for any value 
of the parameter ^. In this case we say that ttq and tti are compatible. Given 
a Poisson pencil tt^ = tti — ^ttq we can often construct a sequence of vector 
fields Yi on M. that have a twofold Hamiltonian form (a so called bi-Hamiltonian 
chain) 



where hi : Ai R are called Hamiltonians of the chain and where i is 




Yi = TTidhi = TTodhi+i 



(2) 
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some discrete index. This sequence of vector fields may or may not truncate 
(depending on existence of Casimir functions). In the case when the Poisson 
pencil TT^ is degenerated but projectable into a symplectic leaf TV (of dimension 
2n) of ttq the bi-Hamiltonian chain on turns into a so called quasi-bi- 
Hamiltonian chain on J\f of the form 

n 

0idH, = eodH,+i + J2 c^v^odH,, i^l,...,n, H^+i = 0, (3) 

where 9i are projections of tt.^ onto J\f, the functions Hj are restrictions of hj 
to J\f: Hj = hj\j\f, and are some multipliers (real functions). And vice 
versa: having a quasi-bi-Hamiltonian chain on the manifold J\f one can 
lift it to a bi-Hamiltonian chain ^ on the extended manifold A4. (Quasi-)bi- 
Hamiltonian chains (called also (quasi-)bi-Hamiltonian systems) possess very 
interesting differential-algebraic properties and are one of key notions in the 
theory of integrable systems, due to the fact that in many cases the systems 
(O and are Liouville integrable T^. Recently, much effort has been spent 
in order to exploit the procedure of solving these systems by the method of 
separation of variables . In this article we will mainly work with the quasi- 

bi-Hamiltonian chains Q rather than bi-Hamiltonian ones, since the pencil 
0^ = 01 — ^00 is always non-degenerated. 

2 Dirac reduction of Poisson bivectors 

We begin by considering the Dirac reduction procedure in a more general setting 
that is usually met in literature. Let tt be a Poisson bivector, in general degen- 
erated on some manifold 7W. Let S = J-'o be a submanifold in a foliation !F of 
the manifold A4 defined by m functionally independent functions (constraints) 
ipi : M —>■ R, i — I, . . . ,m. 

Ts ^ {x e M : ipi{x) = Si, i = 1, . . . , m} 

Thus, iS is a submanifold of codimension m in M. Moreover, let Zi, i — 1, m 
be some vector fields transversal to J-'s, spanning a regular distribution Z in A4 
of constant dimension m (that is a smooth collection of m-dimensional subspaces 
2x C TxA4 at every point x in M). The word 'transversal' means here that no 
vector field Zi is at any point tangent to the submanifold J^s passing through 
this point. Hence, the tangent bundle TM splits into a direct sum 

TM ^ TT®Z 

(which means that at any point a; in we have T^M. = T^^J-s ® with s such 
that X <E Ts) and so does its dual 

T*M = T*T®Z\ 
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where T*!F is the annihilator of Z and Z* is the annihilator of TT. That means 
that if a is a one form in T*J- then a{Zi) = for all i = 1, . . . , to and if (i is 
a one- form in Z* then (3 vanishes on all vector fields tangent to Ts- Moreover, 
we assume that the vector fields Zi which span Z are chosen in such a way that 
dipi, I = 1, TO is a basis in Z* that is dual to the basis Zi of the distribution 

{dipi, Zj) = Zj{(pi) = S^j, (4) 

(this is no restriction since for any distribution Z transversal to T^ we can 
choose its basis so that is satisfied). Finally, let us define m vector fields Xi 
on A4 and functions ipij : ^ M on through 

Xi = Tr{dipi), ipij = {(pi,ipj}j, = {dipi,TTd(pj) = Xj{(pi). (5) 

The functions ifij define an TO-dimensional skew-symmetric matrix ip — [ipij ) , 
i, j = 1, . . .m. It can be easily shown that 

[Xj,X.i\ = Xy^^^^-}^ = TT = TldLpij, (6) 

where {., .I^r is a Poisson bracket defined by our Poisson bivector tt and [X, Y] = 
LxY — X{Y) — Y{X) is the Lie bracket (commutator) of the vector fields X, Y. 

A very special choice of our transversal vector fields Zi originates by taking 
linear combinations of fields Xj with coefficients being the entries of the matrix 

m 

=^((^-i),,X, z = l,...,m. (7) 

Since the constraint functions ipi are functionally independent, the vector fields 
Zi in ^ will indeed be transversal to the foliation T. Moreover, they will auto- 
matically satisfy the orthogonality condition |0J as Zi((pj) = ^Y^=\^'~P~'^')ki^k[.H^ j) 

Let us now consider the following deformation (modification) of the bivector 

tt: 

m 

T^D=T^ -\^X,kZ,, (8) 

i=\ 

where A denotes the wedge product in the algebra of multivectors. This new 
bi- vector tt^i can be properly restricted to Ts for any s G M™ (and thus also to 
S = T{))^ since the image of tt^i considered on a given leaf J-g of the foliation T 
lies in TTs for all s. This is the content of the following theorem. 

Theorem 1 Suppose x Cz J-g- Then for any a £ T*A4 the vector T:]j{a) is 
tangent to Tg i.e. 7r£)(T'*A^) C T^J-g- 

Proof. We will show, that TTuidipk) = 0, k = 1, ...,m, since it means that 
the constraints ipi, i = 1,...,to are Casimirs of ttd (and so are then ipi — Si) 
which obviously implies the thesis of the theorem. Using the definition © of 
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ttd, the obvious fact that (Xi A Zi)dipk = {Xi ® Zi)dipk — {Zi ® Xi)dipk = 
(dipk, Zi) Xi ~ {dipk,Xi) Z^ = Zi{ipk)Xi - Xi{(pk)Zi we have 

TToidipk) = TT{dipk) - ^^Zi{ipk)Xi + \ y^^Xi{ipk)Zi 

i i 

— Xk — ^ ^ SikXi + i ^ ipki{ip^^)jiXj 

= Xk — \Xk — \ ^ SkiXi = 
i 

due to the fact that (fij = —<Pji. ■ 

Theorem 2 The bivector ttd in |3) with Zi as in ^ satisfies the Jacobi iden- 
tity. 

Proof. It is easy to check that our operator ttu defines the following bracket 
on M 

m 

{F,G}^, = {F,G}^ - {F,^,}A^-'hWvG}^, (9) 

(where F, G : Ai ^ M. are two arbitrary functions on A4) which is just the well 
known Dirac deformation (IHj of the bracket {., .^t^ associated with tt, and as 
was shown by Dirac fj] it satisfies the Jacobi identity. ■ 

Remark 3 IfC:A4^Misa Casimir function ofir, then it is also a Casimir 
function ofno, since in this case yields 

m 

{F, C}^^ = {F, C}^ - J2 ^^}AV~%WJ,C}^ = 0-0 = 0. 

We also know from Theorem 01 that the constraints (pi are Casimirs of the de- 
formed operator ttjj. Thus, we can informally state that Dirac deformation 
preserves all the old Casimir functions and introduces new Casimirs ipi. 

It is now possible to restrict our Poisson operator ttd (or our Poisson bracket 
{., .Itto) to a Poisson operator (bracket) on the submanifold S (i.e. a symplectic 
leaf ifi — ... — ifm = of TTo) in a standard way (reduction to a symplectic 
leaf). Namely, for arbitrary functions f,g : S R one defines the reduced 
Dirac bracket {/, g}i^ on S (with the corresponding Poisson operator tt^) as the 
restriction of the bracket of arbitrary prolongations of / and g to A4, i.e. 

{f,9}R = {F,G}D\s 

where F,G : A4 ^ M. are two functions such that f = F\s and g = G\s. 
Of course, an identical construction can be induced on any leaf JF, from the 
foliation J-. There arises, of course, a question how 'robust' this construction 
is, i.e. to what extent the obtained deformation no and the reduction ttji are 
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independent of the choice of particular functions that define our submanifold 
S. This issue will be partially addressed in next sections. 

The results of this section as well as the results presented in Sections 4 and 
5 suggest that the concept of the classification of constraints as being either 
of "first-class" or of "second-class", proposed by Dirac, should be reexamined 
when one looks on the problem from the point of view of Poisson geometry. 
First of all it is clear that the procedure called Dirac reduction has two different 
levels. The first level we call Dirac deformation as we deform a Poisson bivector 
TT from manifold M to another Poisson bivector ttd on the same manifold A4. 
A sufficient condition for the existence of tt/j for a given set of constraints ipi, 
i = 1, ...,m is a nondegeneracy of the Gram matrix if. From the construction 
all constraints tpi are Casimirs of Tr^j.The second level of the construction we 
call Dirac restriction as we restrict a Poisson bivector ttd to its symplectic leaf 
J^s- In the original Dirac construction it was a particular one, namely J^o = <5. 
A Poisson bivector on is denoted by ttr. For the existence of the second 
step additional restrictions on ipi have to be imposed. Actually, tth has to 
be nonsingular. In a standard classification it means that we have to exclude 
constraints of first-class. Let us remind that a constraint (pk is of first class if 
its Poisson bracket with all the remaining constants ipi vanishes on S, that is if 



Otherwise ipk is of second-class. In general, first-class constraints make tt^ sin- 
gular so that the Dirac reduction procedure can not be performed. Nevertheless, 
the condition (|10|l seems to be too strong. In sections 4 and 5 we demonstrate 
situations where constraints are of first-class but the singularity is 'removable' 
and so the Poisson bivector tt/j is well defined. It suggests that the classification 
of constraints given by Dirac should be reformulated in the context of Poisson 
pencils and Poisson geometry. 

3 Separable quasi-bi-Hamiltonian chains of Be- 
nenti type 

In the following section we briefly remind basic facts about separable Hamilto- 
nian systems on Riemannian manifolds, which form a special class of quasi- 
bi-Hamiltonian chains 0,^3 1 known also as the so called Benenti systems 
|23,|22|- Let {Q,g) be a Riemannian manifold with covariant metric tensor 
9 = (ffy) ^'^d with the inverse (contravariant metric tensor) = G = (G"-'). 
Let {q^, g") be some coordinate system on Q and let {q^ , ...,q^,pi, ...,Pn) be 
the corresponding canonical coordinates on the phase space Af = T*Q with the 
associated Poisson tensor 



where J„ is n x n unit matrix and 0„ is the n x n matrix with all entries equal 
to zero. Let us consider the Hamiltonian E : M ^ M. ior the geodesic flow on 



{(Pk,fi}-^\s = 0, i = l,...,m. 



(10) 




(11) 
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Q: 

n 

As it is known, a (1, l)-type tensor B = (_Bj) (or a (2, 0)-type tensor B = 
{B^^)) is called a Killing tensor with respect to 5 if { J^i^Cy^ piPj , i?}^ = 

(or {J^i^T'' PiPj J — 0)- important generalization of this notion is 

formulated in the following definition. 

Definition 4 Let L = (L*) 6e a second order mixed type (i.e. (l,l)-j tensor 

n 

on Q and let L : Af -^M. be a function on M defined as L — ^ {LCy-'piPj , 

n 

where LG is a (1, 1) tensor with components {LGY^ = L\.g^^ . If 

{L,E}e„^aE, where a=f]G^^|4p„ f^Tr{L), 

then L is called a special conformal Killing tensor with the associated potential 
f = Tr{L) ini. 

The importance of this notion lies in the fact that on manifolds with tensor 
L the geodesic flows are separable. There exist, in this case, n constants of 
motion, quadratic in momenta, of the form 

n n 

Er=Y^ Al^p^p, = {KrCrp^pj, r = 1, n, (13) 

n 

fe=i 

where and if,- are Killing tensors of type (2, 0) and (1. 1), respectively. More- 
over, all the Killing tensors Kr are given by the following 'cofactor' formula 

n-l 

co/(e/-L) = ^X„-,e, (14) 
i=i 

where cof{A) stands for the matrix of cofactors, so that cof{A)A = {deiA)I. 
Notice that Ki ~ /, hence Ai = G and Ei = E. Since the tensors Kr are Killing, 
with a common set of eigcnf unctions, the functions Er satisfy {Es,Er}Bg — 
and thus they constitute a system of n constants of motion in involution with 
respect to the Poisson structure Oq. So, for a given metric tensor g, the existence 
of a special conformal Killing tensor L is a sufficient condition for the geodesic 
flow on A/" to be a Liouville integrable Hamiltonian system. 
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The special conformal Killing tensor L can be lifted from Q to a (1, l)-type 
tensor on A/" = T*Q where it takes the form 

"-{fI-). n--l-S^v),~-^^p'L),. (15) 

The lifted (1,1) tensor N is called a recursion operator. An important property 
of N is that when it acts on the canonical Poisson tensor it produces another 
Poisson tensor 

ei = = _° T 

compatible with the canonical one (actually 6q is compatible with N^Oq for any 
integer k). 

It is now possible to show that the geodesic Hamiltonians E^- satisfy on 
M — T*Q the set of relations 

9odEr+i ^ 9idEr + PrOodEi, En+i^O, r = l,...,n, (16) 

where the functions Pr{q) are coefficients of the characteristic polynomial of 
L (i.e. minimal polynomial of N), which is a special case of the quasi-bi- 
Hamiltonian chain Q pp. 

It turns out that with the tensor L we can (generically) associate a coordinate 
system on JV in which the flows associated with all the functions Er separate. 
Namely, let {X^{q), A"(g)) be n distinct, functionally independent eigenvalues 
of L, i.e. solutions of the characteristic equation det(^J — L) = 0. Solving these 
relations with respect to q we get the transformation A — > g 

q'^a,{X), i=l,...,n. (17) 

The remaining part of the transformation to the separation coordinates can 
be obtained as a canonical transformation reconstructed from the generating 
function W{p,X) = ^iPiCtiiX) in the standard way by solving the implicit 
relations pi — ^^^Z'^'' with respect to pi obtaining pi = (3i{\,p,). In the (A,/i) 
coordinates, known as the Darboux-Nijenhuis coordinates (DN), the tensor L 
is diagonal L = diag(A"'^, A") = A„, while the Hamiltonians H13|) of our quasi- 
bi-Hamiltonian chain attain the form [3| 

^^(A,M) = -1^^^— r = l,...,n, 

where 

A, - n (^'' - ^')' 

k—l,...n, k^i 

Pr{X) are symmetric polynomials (Viete polynomials) defined by the relation 

71 

det(a - A) = (^ - X'm - A')-(C - A") = J2 p^^^^ (18) 

r=0 
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and where fi are arbitrary smooth functions of one real argument. 

It turns out that there exists a sequence of generic separable potentials Vr''^ , 
A: e Z, which can be added to geodesic Hamiltonians Er such that the new 
Hamiltonians 

Hriq,p)^Eriq,p) + V}''\q), r = l,...,n, (19) 

are still separable in the same coordinates (A,/i). These generic potentials are 
given by some recursion relations 01,501 ■ The Hamiltonians Hr : TV ^ R in 
p9(l satisfy the following quasi-bi-Hamiltonian chain 

9odHr+i ^ OidHr + PreodHi, i?„+i = 0, r = l,..,n. (20) 

In the DN coordinates the Hamiltonians Hr attain the form 0] 

Hr{\ f^) = - 2^ -gj- ^- , r = l,...,n, (21) 

i—l * 

where potentials Vr'''' enter Hr as 7i(A*) — (A*)""'"'^"^. From it immediately 
follows that in (A, n) variables the contravariant metric tensor G and all the 
Killing tensors Kr are diagonal 

C - , (Kr)^- g^^d^. 

Moreover, in the (A, /i) coordinates the recursion operator and the tensor 9i 
attain the form 

/ A„ 0„ \ f On K \ 

\ 0„ A„ \ -A„ 0„ ) 

while ^0 remains in the form Hll() since the transformation ((?,p) — > (A, /x) is 
canonical. 

The quasi-bi-Hamiltonian chain H2U|) on M can easily be lifted to a bi- 
Hamitonian chain Q on the extended manifold M. = T*Q x R = Af x II 

m 

Having such a complete picture of separable quasi-bi-Hamiltonian chains on 
Riemannian manifolds, one can ask a question: what kind of holonomic con- 
straints can be imposed on considered systems so that this quasi-bi-Hamiltonian 
separability schema is preserved? The simplest admissible case of such con- 
straints will be considered in the next sections. 



4 Reduction A" = 

Let us consider a particle moving in our Riemannian manifold Q equipped with 
the coordinates (g^, ...,g"). Suppose that this particle is subordinated to some 
holonomic constraints on Q defined by the set of relations 

Mq)=0, k^l,...,r (22) 
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that define some submanifold of Q. The velocity v — X^i"^*^ '^^ ^^lis particle 
must then remain tangent to this submanifold so that 

= (dipk,v} ^J2-£r'''- 

and thus in our coordinates = G^^pj the motion of the particle in the 
phase space M = T*Q is constrained not only by the r relations but also 
by the r relations 

^r+k{q,p) ^ J2 G'^{q)^^p,^0, k=l,...,r. (23) 

that are nothing else than the lift of H22I) to J\f. We will call the constraints 
l|23(l a g- consequence of the constraints H22(l . as they are natural differential 
consequences of (I22II at given metric tensor g. The constraints (|22|I - H23(I define 
a submanifold 5 of A/" of dimension n — 2r. 

Let us now consider our quasi-bi-Hamiltonian chain H2U|) in J\f. We would 
like to know what types of holonomic constraints (|22|l on Q do not destroy the 
separability of the constrained chain? This is a complicated question and due 
to its nature it is most convenient (for a moment) to consider it directly in 
our separation coordinates (A,/x). Thus, in this section we will analyze only 
a very special choice of the functions ipk in (|22(l . Namely, we put r = 1 (so 
that n — 2r — n — 2 which corresponds to to = 2 in Section 2) and define the 
corresponding function ip[ in (A, /i) variables as 

<p'i(A) = A" (24) 

Since the metric tensor G — (G*^ ) has in (A,/i) coordinates the diagonal form 
G = diag (/i(A*)/Ai) , i = 1, . . . , n and since the equation 1)23(1 has an invariant 
form (i.e. in (A, /i) coordinates it has the same form with q and p replaced by A 
and respectively) the g-consequence of reads as 

^2(A,M) = ^^^Mn (25) 

(we use ' here since soon we will modify the constraints (|24|I - H25() to a simpler 
form). These two constraints define a subset S' = {f'l = 0,(P2 — 0} of TV. We 
will however restrict us to a submanifold S = {A" = 0, /i„ = 0} neglecting the 
term /„(A")/A„ in (for 7^ 5 C 5' and in case of the systems when /„ 
never vanishes S and 5' coincide) . The reason for this is that various calculations 
with the help of these new constraints A" = 0,/i„ = are simpler. Let us now 
perform the Dirac reduction of our Poisson operators 6q and 9i 



70 = 



0„ In \ n ^( 0„ A„ 

-/„ 0„ j ' ^ I -A„ 0„ 



(26) 
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and of the corresponding quasi-bi-Hamiltonian chain (geodesic case) or 
1)2011 (potential case) on the submanifold S. We can do it either by using the 
constraints (|24|I - H25|) or by the constraints that directly describe S: 



^i(A) = A" = , ^2 (A, A*) 







(27) 



(of course ipi = ip'i). It turns out that the corresponding Dirac deformations 
(given by ^ or equivalently by lO) do,D and 9'q ^ of 6*0 are not equal. Similarly, 
Qi.D 7^ ^'i D- However, after reducing Sq.u and 6'^ jj on S we obtain the same 
Poisson operator (and similarly with 9i_d and 9'^ jj). 



Theorem 5 With the notation as above, the reduced operators 9i ji and i 
S satisfy 9q = 9q ji and 9[ — 0i,_r. // we parametrize S by the variables 



i,R 



(A\. 



,A"- 



70, i? 



. , ^„_i) then 

0,1-1 In-l 
— In-l 0„_i 



On-1 

-A„_i 



A„-i 
0„-i 



(28) 



with A„_i = diag(A^, . 



,A"- 



Proof. Let us start by calculating 0o,d- In this case the 2x2 Gram matrix 
f — iWi, ^j}0g) is of the form 





-1 



and thus according to (0 we have Zi = X2 and Z2 = —Xi. The formula JHJ 
yields 

d d d _ d 



70,15 



9q — Xi A X2 



aA" 



^ 9A» 

i—i 



since Xi = 9odipx = —djd^in and X2 = 9odip2 = d/dX 
coordinates we have 



70,13 



o„ 


In-1 0(„_i)xl 
Olx(n-l) 


Oix(„-i) 





. So, in (A,/i) 



(29) 



The operator l|29|l has thus the n-th and the last column and the n-ih and the 
last row filled with zeros. These two rows and two columns of zeros correspond 
to the fact (see Theorem ^ that the constraints ipi = A" and (^2 = are 
now Casimir functions for 9q d. We can thus directly project 9o £i onto the 
symplectic leaf S = {A" = /in = 0} of ^o.d by simply removing the two zero 
columns and zero rows from (|29|) . which yield the operator 9o^r that written in 



coordinates (A"'^, 



• Mn-i) Oil '5 attains the form as in H28|l . Similar 



computations show that in the case of 9i j) 



" A" ■ 









-A" 
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and thus (see CO))) (^i and ip2 are first-class constraints. Indeed, following (Q 
we have Zi = pr^2 and Z2 = —^Xi and the formula (jSJ yields 

Oi,D = di — l^^i ^ ^2 

which seems to be singular on S. Nevertheless this singularity is 'removable' 
since Xi = Oidipi — — A"9/9/i„, X2 — Oidif2 ~ A"9/9A" and hence 



■t— 1 



9A" ^ dX 

l—l 

So, the matrix form of 9i^d in (A,/i) coordinates becomes 



o„ 


A„-i 0(„_i)xi 

Olx(n-l) 


-A„_i 0(„_i)xi 
Oix(„-i) 


o„ 



and the projection of this operator onto the symplectic leaf S = {A" = /i„ = 0} 
of yields exactly 0i^fi as in if^ . 
Passing to (p'l , (/92 notice that 



¥^'12 



/«(A") 



so one can have a situation when <y3'i2|5 = which again leads to first-class 
constraints according to Dirac classification. Nevertheless, calculations similar 
to these for ipi, ip2 show that singularities are again 'removable' and we end up 
with the same form of Oq d and 9i^d- ■ 

Of course, the reduced operators 0o,-R, and 61^11 are compatible. Moreover, it 
is easy to see that the Hamiltonians ()21|l restricted to S become 



Hr\s=-Y: 



9A» 



r = l,. 







where /,,fl(AO = /*(A0/A' , -/^A>^') = 7.(A0/A\ z = 1, . . . , n - 1 and pr^W 
are Viete polynomials of dimension n — 1. Having all this is mind, we can 
formulate the following corollary. 



Corollary 6 Both Dirac reductions ^^-^3) and \21^ of the quasi-bi-Hamiltonian 
system \2(]\) onto S lead to the same quasi-bi-Hamiltonian system of the form 

do,RdHr+i^R = 9i^RdHr^R + Pr,Rdo,RdHi^R, HnM — O, r = l,...,n— 1. 

which is thus separable in the variables (A^, . . . A"~^, /ii, . . . Pn-i) that coincide 
on S with 'first' n — 2 separation coordinates from our coordinate system (A,/i) 
that separates \2(]\) . 
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Because the reduced quasi-bi-Hamiltonian chain has exactly the form (|2U|I . 
it can be put into a bi-Hamiltonian chain on 5xR. 

Of course this procedure of Dirac reduction, in principle could be performed 
directly in {q,p) coordinates, since the formula (|SJ) has a tensor character and 
thus yields the same result no matter what coordinate system we choose for 
actual calculating of deformations Oq^d and Oi^u. There is however a major 
obstacle here: we are usually not able to express the constraint A" = in 
'physical' coordinates {q,p), as equations (|17|) are noninvertible in general, i.e. 
there is no algebraic way of solving them with respect to A (if we could, then the 
second constraint could be computed as the g-consequence of the first one 
calculated directly in {q,p) coordinates, with the help of Thus, although 

the picture presented in this section is clear, nevertheless it is somehow useless 
once we are given the quasi-bi-Hamiltonian chain in natural {q,p) coordinates. 
In the next section we will demonstrate how this problem can be defused by the 
use of the so called Hankel-Frobenius coordinates. 



5 Reduction in Hankel-Frobenius coordinates 



As it was demonstrated in the previous section, the constraints ipi = A", (p2 = 
jjLn preserve the separability on the constrained submanifold S but are very 
inconvenient to handle with, as in general we do not know how to express them 
in the original coordinates {q,p). More convenient for this purpose is the set 
of so called Hankel-Frobenius coordinates {pi,Vi)^^^ They are non- 

canonical coordinates, related to the separated coordinates in the following way 



where 



dX'' Afc 
(Ai)"-i ••• Ai 1 



1, ...,n. 



(30) 



k 



Vn 



(A«)«-i ... A" 1 

is a Vandermonde matrix and pi are Viete polynomials (|18|l . The mapping H30|l 
is not a point transformation on Q and therefore it makes no sense to distinguish 
covariant and contravariant indices now. Applying the map H3U|I we obtain the 
form of operators 6*0 and 6i and the 2ri-dimensional recursion operator iV„ (|15|l 
in (p, v) coordinates: 

(Q ... 1 \ 
... 1 pi 



='0 = 







Un 





V 1 pi 





-FnU„ 



F U 




F U 



Pn^l I 

- J] f'^ 



(31) 



(32) 
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Fn 

F„ 



-Pi 

-P2 



\ 





-Pn-l 1 

V -Pn / 

We will now consider yet another deformation ||SJ| of Oq and 6i, given in (A,/i) 
variables by the constraints 



P„(A) = (-1)"A1a2...A" 
= «„(A) = {-ly^Y^tLx^X^ . . . A^-iA^+i • • • A" 

.7 = 1 ^' 



(33) 



equal to last pair of Hankel-Frobenius coordinates. The constraints (|33|l are 
related with the constraints lf77|) as 



(34) 



where the functions 

V'2 



(-i)"aia2---a"- 



+1 



V n — 1 



and ■03 = (— 1)"A"'^A^ • • • A"^^/A„ never vanish at 5 = {(pi — ^2 — 0} (in fact, 
4'2\s — ~1 and ^p3\s — —{p-i + ■ ■ ■ + ) and thus the constraints H33(l define 

(locally) the same submanifold S as the constraints (I27|l . The corresponding 
deformations (?q ^ and 9'{ of 0o and di will of course not be equal to Oq^d and 
0i^D, but again it turns out that their reductions on S will coincide with the 
corresponding reductions of Oq^d and Bi d on S. 



Theorem 7 In the notation as above, 9q ^ — Oq_r and 0q — Oi 



OM 



Proof. We will use rather than ((SJ since it turns out that the calculations 
are in this case simpler when one uses bracket definition of Dirac deformation 
than bivector definition. Applying Q, we easily get that for any two functions 
A,B:Af ^ R 



{A 'P2}e„ {B, ipi}g^ - {A, ipi}g^ {B, ip2}g^ 



{<Pi></'2}eo 

where of course {(fi, f2}g^^ = 1 and so it does not vanish on 5. Similarly 

{A, ^'i). {B, ^'i). - {A, VP'/}. {B, ip'i). 



{A,B}, 



(35) 
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Using the relations H34|l between the deformed constraints tp'l and the original 
constraints (pi, the Leibniz property of Poisson brackets and the fact that ipi 
and ipz depend only on A, we obtain 

+ l/'lV'3{V'l,'^2}0„ + V3i¥32 {V'l,'03}eo, 

so that {ip'l, ip2}eg \s = '^iV'a {(y^"!, V2}0(, Is = '^li'sls- Similar calculations show 
that 

and thus the factors ipi'4'3 i^i the numerator and in the denominator of H35|l 
cancel and we conclude that 



which is the same as to claim that 9q 



The proof that 0" i? = ^i,-R is 



similar: first one shows that 

{fl:V2}ei Is = 1pl1p3Wl,f2}g^ Is = V'l'/'slsA" 

which is by the way equal to zero on S. However , we also get 
{{A, ^'i},^ {B, ^'{},^ - {A, ^'{},^ {B, ^^'}, J 1^ 
= i/'iV'a {{A, ip2}e, {B, tpi}g^ - {A, ipi}g^ {B, (^2}eJ | 



so that on S 



V'lV's {{A, ip2}e, {B, tpi}g^ - {A, ipi}g^ {B, tp2}gj 



{Ay2}g^ {B,tpi}g^ - {A,tpi}g^ {B,tp2}g^ 



A" 



since the term A" in the last expression does not cause any singularity: for 
every possible combination A, B — , the numerator in the above expression 
is either zero or some multiple of A". This proves that 6'^ ~ Oi,r- ■ 

Before we proceed with the main theme of this article, let us make a digres- 
sion: we will establish the form of Oq^r and ^i.jj in (p, u)-coordinates. 



Lemma 8 In (p, v)- coordinates we have 



%.R - {Nn-lY 



C/„-i 

-Un-i 



0„-l 



n,R 



CX„_i 

-Un-l 
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This lemma can be proved either by direct calculating of the deformation 
formula ^ in (p, ti)-coordinates or by transforming both Oq^r and Bi^r with 
the map that is the restriction of the map (|30|) to S. This restricted map has 
no longer the form H3U|) and this is why the operators 6*0, ^ and 6i^r does not 
transform respectively to the operators of the form and H32|) with n replaced 
by n — 1 . 

Having established Theorem[71we can now, for a given (quasi)-bi-Hamiltonian 
system (|2U|I , perform the Dirac reduction A" — fin = directly in physical coor- 
dinates {q,p) by performing the equivalent reduction pn = f « = provided that 
we are able to express the constraints (I33II directly in physical coordinates. The 
first of the constraints in i.e. p„ = is given directly in {q,p) coordinates 
once the system (|16|) has been given (this is the reason for which we consider 
this constraint and Hankel-Frobenius coordinates) . Below we shall show how to 
express the second constraint {vn = 0) in physical coordinates. 

Let us first observe that in case of systems for which fi{X^) — const (that 
does not depend on i) the constraint w„ = calculated in (A, fj.) coordinates is 
just the (/-consequence of the first one (p„ = 0) , so that in this particular case 
we can easily get the function 4>2{q,p) by calculating the expression with 
an appropriate metric tensor g that can be found in hi. When the functions fi 
are more complicated, we must proceed differently. 

Theorem 9 Assume that /i(A*) — /(A*), i.e. that the functions fi do not 
depend on i. Let also F he the matrix given by Then the coordinates Vi 

expressed in {q,p) coordinates have the following form 




with F being expressed in q. 

Proof. Let us expand the function / in a formal Laurent series 

/(A)=^a„,A", a„eM , m e Z/ C Z. (37) 

■m 

Let also , A: = 1, . . . , n be the ^-consequence of pk so that 

<(A,A*) = E ^''^ '^^^'^'^^ = G^'il)fjPs- (38) 

i—l * — 1 

Setting n = (^i, . . . , PnV , v = (wi, . . . ,w„)^, v' = {v[, . . . ,w^)^ and using the 
fact H3()|l that /j, — Vv (where V is the Vandermonde matrix) we easily obtain 
that 

i;' = y-Miag(/(Ai),...,/(A"))y« 

Moreover, from the transformation between (A, fi) and (p, v) coordinates one 
finds that 
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so that due to H37|l 



^' = l/-i^a„ diag((Air,...,(A"r) Vv^V~'[Y,arn A'" Vv 



This relation is valid in every coordinate system and thus v{q, p) — f{F{q))v' {q, p) 



In case when the assumption of the theorem is not satisfied, we can not 
express w„ in {q,p) variables using the method presented in the proof. Since in 
practice this situation is very rare, we choose not to discuss it in this article. 
Notice that if / = const then f{F) = const I and the formula H3t)|) for reduces 
to the ^-consequence of p„, as it should be. 

Remark 10 One can ask the question why, in the case of f ^ const, instead of 
taking = Vn o,s above, we do not choose if 2 to he simply the g- consequence of 
= Pn, o,s both pairs of constraints describe the same S. The answer is due to 
one of the fundamental observations of presented paper, mentioned at the end 
of section 2. Actually, a constrained submanifold S can be defined by infinitely 
many different pairs of constraints, all of them of the first class according to 
the classical Dirac classification. That is, either the Gram matrix is singular: 
{if I, }e ~ {'^1 ' ¥'2 }e Is = 0, for one or both 9i, hence, in principle, there 
is no Dirac deformation 9i_D or the restriction of 6i,D on S is not possible. 
However, among all these pairs of constraints, there are exceptional pairs, like 
(A",/i„) or [prnVn), for which singularities are 'removable' and Dirac restricted 
Poisson pencil is nonsingular. So in fact, these particular constraints have to 
be considered as second order constraints. 

Notice also that in order to determine the constraint u„ from (|36|l we have 
first to find the function / for example from Ei, written in (A, /i) coordinates. 
So in practice, we cannot avoid the calculation of separation coordinates, and 
this is the price we have to pay. Nevertheless, in general, it is not difficult to 
find the transformation (|17|l and then the geodesic Hamiltonian Ei{X, fi) for the 
original system. Also, one has to bear in mind that even though the reduction 
procedure looks trivial in (A, /i)-coordinates, it is usually not possible to express 
the obtained reduced system back in ((7,p)-coordinates since our Hamiltonians 
Er^R have different function / than the Hamiltonians Er. That is why we had 
to use Hankel-Frobenius coordinates. 

6 Examples 

In this chapter we will illustrate the introduced ideas by two examples. In 
our first example of Dirac reductions of separable systems we will consider the 
so called first Newton representation of the seventh-order stationary flow of 



m 




which due to yields (jSHl- ■ 
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the KdV hierarchy [201 > 121 ■ It is a Lagrangian system of second order Newton 
equations 

10(qi)2 + V 

20q^q^ - 8(g2)2 + 30{q^)^q'^ - 15(9^)"* + c 

(where the subscript , t denotes the differentiation with respect to the evolu- 
tion parameter i G M). The above system can be represented as a quasi-bi- 
Hamihonian system on TV belonging to a quasi-bi-Hamiltonian chain of the 
form H2U|) with n = 3, with Hamiltonians 



9,tt 



Hi 

H2 



PlP3 



\q^pl - WpI 

h 



10(g 



q-" - Aq^q^ + ^q^{q^f - lQ{q^fq^ + iiq^f 



\(f'P2Pz 



^\{q'f-2{q'f + {q^f 
H,^\{q'fpl + \{q'fpl + 



2P1P2 
1^2 „3 



WpiP,+2iq'nqr + W)v 



Gq'q^q 



(39) 
(40) 



\q^piP2 



jq^pips 



112 

4? g P2P3 



\q^P2P3 



i{q'f{q^f+qHq^f + Wfq^ + 2q\q^f 



+ |(9^)V + (gW-(9^)V9^ 
with the corresponding canonical operator 



((TT|l and 61 of the form 









-1 











q^ 


-1 














-q^ -2q^ 


P2 


P3 


1 





-P2 








1 -q' 









+ P3 = 




4^ ■ 



From the form of Hi one can directly see that the inverse metric tensor G 
expressed in (g, p) variables has in this example an anti-diagonal form 



1 

1 

1 



(41) 



while the conformal Killing tensor L has the form 

^1 -1 



2 



q 

q' 
2q3 



-1 
2 1 

q q 



which substituted in H14|l yields the geodesic parts of all the Hamiltonians (|4()(l . 
Our quasi-bi-Hamiltonian system turns out to be separable in (A, /i) coordinates 
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defined as above and it turns out that in this case fi — 1/8 = const and 
7i(A') = 16(A*)^ in (|21() . Thus, we can easily find out the constraints ps = 
V3 — directly in {q,p) coordinates. From the form of H3 in (|40|) one can 
see that (p'( = Psiq) = —j{qiq2 + 93) and that = v^{q,p) is just the g- 
consequence of p^. An easy computation of with the use of H41|) yields that 
i'3('?:P) = ~\{Pi + P2q^ + P3q^)- Performing the Dirac deformation ^ of 9o 
and 9i with the use of constraints ip'( and 1^92 yields 



yo,D 



and 











q' + iq'r 


-q' 




-1 










1 9 

~q q 






-q' 


1 








-{q'r ' 


-q'q' 


q^ + 


{q'f 


2q^ + (gi)2 









q^P2 


- q^P3 P2 








q'^pa - q^p2 







P3 








-P2 




-P3 
















-1 















g2 





-1 












-2q^q^ 


q' 


q' 




ei,D - 1 






2q^q^ 


P2 


P3 






1 





-q^ -P2 















1 


-q^ -P3 











so that Oi^D differs from 61 only at entries (3,4) and (4,3), where q'^ was deformed 
to —q^q^ (notice that on S indeed we have q'^ = —q^q'^). We will now pass to 
the Casimir variables chosen as variables 



(9\ 9^ , </'i (9) , ¥'2 (?> -P) ; -P2 , P3 ) 



(42) 



since due to the fact that it is easiest to eliminate q'^ and pi from the system of 
equations = Lp'l{q), ^2 = v'^iq^P) ^i^l parametrize our submanifold by 
the coordinates {q^,q^,p2,P3). In the variables H42|l the operators Oi^o attain 
the form 



yo.D = 



2g2 + (gl)2 



and 



n.D 



















-q' 


-1 


















2g2 


-q' 
















































q' 


2q' 













P3 




1 


q' 










-P3 



















-1 
























-1 
















































1 
























1 



















respectively. Observe that in these new variables all the entries at third and 
fourth rows and columns are zero, as it should be, since now the Casimirs ifl 
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and are part of our coordinate system. We can now write down the reduced 
operators 0o,fl and Oi^r in variables {q^,q^,P2,P3) by simply removing the zero 
rows and columns (we would also have to put = but our matrices Oi^D do 
not contain any variables (p'^ in their entries). This yields 



yo,R — 











-q' 


-1 


1 








2q' 


-q' 


2g2 + (gl)2 


q' 







P3 




1 


q' 


-P3 






yi,R 



O2 



h 
O2 



so that Oi^R attains in the variables (q^,q^,p2,P3) the canonical form. An easy 
calculation yields that the Hamiltonians Hi restricted to S become 

Hi,R{q\q^P2,P3) ^Hi\s = -q'p2P3 + \pI - q^pl - '20{q^)^q' 
+ 12gi(g2)2 + 3(^1)5 

H2Mq\q\P2,P3) = H2\s = {q^ + P2P3 + %q^f{q^f 



-H"3,i^(g^g^^'2,P3) 



2 ^ 







while the functions pi^R are: pi^R = —q^, P2,R. = j{q^)'^ + ^q^, P3,R ~ 0. Thus, 
we have obtained on 5 a new separable quasi-bi-Hamiltonian system of the form 
l|2()(l with n = 2. This concludes our first example. 

As a second example we will consider a separable (quasi)-bi-Hamiltonian 
system with a non-trivial function /. More specifically, we shall consider a 
Lagrangian system of second order Newton equations 



q% 



6V +48((7^)2 +4C + 24 



16q\l+q^) 
q%^m'+A{q'f 
q3, = 4q3(3 + 4g2) 

This system is a part of a quasi-bi-Hamiltonian chain with Hamiltonians 

H^^pI+pI+pI~ 4(ql)2(l + q2) _ (^3)2 (3 ^ 4^2) 

-4(3<z2+4(g2)2 + 2(g2)3) 
H2 = - (2g2 + 1)pI + 2gViP2 -pI + 2q^p2P3 - 2pl{q^ + 1) 

- {q^r {{q^ + 2iq'r + ^q'f + V + 2) - {q^ 
-f-8(g2)3 + l6(g2)2 + i2g2 

H3 = -(g')'p? + 2q'q^p,p3 - 2q^p2P3 + (l + 2q^ - {q'f) pi 
+ {iq'f+'iiqY + iq'r + &q' + S) {q')', 



A{q 



2\2 



(43) 
2q^) 
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with the canonical Poisson tensors 9o and 












1 


q' 














q' 


2q^ + l 


q' 














q' 





-1 










-Pl 







-2g2 _ 1 


-q' 


Pi 





P3 





-q' 








-P3 






and with the functions 

Pl{q) = -2{q^ + 1) , P2{q) = 1 - - iq'f + 2q' , pM = [q^f 

The above system is a particular example of so called bi-cofactor systems |27| . 
|28| . \2'^\ . |29j . pTHj separated recently in il^j. In Darboux-Nijenhuis coordinates 
(A,^) it turns out that now /(A*) = 4A*(A* - 1) and 7i(A') = (A')^ This means 
that in order to find the constraint (p2{q,p) we have to use the formula (|36|l with 
fiiq) = (9^)^- Plugging these functions into H3t)|l and using the above forms of 
G and F we get a rather complicated expression for 



'P2{q,P) = V3{q,p) 



Pl- 



3^2 



P2+P3 



2g2 + (gl)2_l 



Notice that the constraint Lp'^ seems to have a singularity on the surface = 0. 
However, it turns out that on the submanifold S = {(fi = <^2 = 0} we have = 
and the singularity disappears (see below). Having obtained the constraints 
ip'l and in [q^p) variables we now proceed as in the first example. First we 
find 



JO,D — 



(where Q are some rational functions of q and p that do not vanish on S) and 












1 





Ci 














1 


C2 




















-1 














Ca 





-1 













-Ci 


-C2 





-C3 


-C4 









































-1 


-q' 







q' 


-2g2 - 1 










-q' 






Pi 



2g2_ 




Pi 1 



ill 





-2P3 





q' 




2P3 





(notice that in both cases is a Casimir function). It can be checked that the 
Schouten bracket of 0o,d and 9i^d is equal to zero: [Oq^d , ^i,-d]s = so that 
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the operators 0o,d and 9i^d are compatible. Before projecting 6'o,d and 9i d 
onto S we perform a transformation of variables to Casimir variables 



{q^,q^, fi (q) ,pi,P2, (P2 {q,p)) 



(44) 



since this time it is easier to eliminate and from the system of equations 
(p'( = ipi{q), = ^2{q-<p)- Explicitly we have 



q 

P3 



(45) 



Vi'- 



q^P2 



2q^q^ 



iq') 



l^3 



As we have mentioned, we have now P3I5 = 0. We will thus parametrize our 
submanifold S by the coordinates (g^, q^,pi,p2). After expressing the operators 
9i^D in the variables H44() with the help of relations (|45() we can easily project 
these operators on S. As the result we get a new (reduced) quasi-bi-Hamiltonian 
chain of the form with n = 2, which is the variables {q^ ,q^ ,pi,p2) is 
determined by 



O2 h 
h O2 






-1 



q 


-2g^-l pi 



q 

2g2 + l 

~Pi 




and by the restricted Hamiltonians Hi 



H,.n = H,\s=pI+pI- 4(gl)2(l + q^) - i {3q^ + 4{qr + ^iq^) 



H. 



2,R 



-(gi)2((gi)2 + 2(g3)2+4(g2)2 



2 
P2 



H: 



3.n. 



+s{qy 

H3\s = 



+ 16(g2)2 + I2q^ 



while the functions pi^R are given by pi^R = —2c{q^ + 1), p2,R. = 1 + 2q^ — (g^)^, 
Pz,R = 0. 

Notice that when we take (p2{q,p) = 2q^p^, i.e. the differential consequence 
of ip'{{q,p) = {q'^y, then {t/c" , </'2 l^i = and the Dirac deformation 9i^d does 
not exists (see Remark [Tn|) . 



7 Conclusions 

In this article we have focused on the problem of Dirac deformation and Dirac 
reduction of Poisson operators and Poisson pencils. We have presented a proce- 
dure of performing Dirac reduction in quasi-bi-Hamiltonian systems of Benenti 
type that do not destroy the separability of these systems and that moreover 
do not change the separation variables. The method presented is not, however, 
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general in the sense that it does not provide us with ah such reductions, but 
only with their subclass. Moreover, for the moment the procedure works only 
for those systems, for which the functions fi are all equal, i.e. when /; does not 
depend on i. Last but not least, the presented procedure is valid only in case 
when one of the Poisson structures of the system is canonical - the non-canonical 
case must be studied separately. These issues will be addressed in a separate 
paper. 
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